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Abstract 

d- and /-shells support a large number of local degrees of freedom: dipoles, 
quadrupoles, octupoles, hexadecapoles, etc. Usually, the ordering of any mul- 
tipole component leaves the system sufficiently symmetrical to allow a second 
symmetry breaking transition. Assuming that a second continuous phase transi- 
tion occurs, we classify the possibilities. We construct the symmetry group of the 
first ordered phase, and then re-classify the order parameters in the new symme- 
try. While this is straightforward for dipole or quadrupole order, it is less familiar 
for octupole order. 

We give a group theoretical analysis, and some illustrative mean field cal- 
culations, for the hypothetical case when a second ordering transition modifies 
the primary T xyz octupolar ordering in a tetragonal system like URu2Si2. If 
quadrupoles appear in the second phase transition, they must be accompanied 
by a time-reversal-odd multipole as an induced order parameter. For O xy , O zx , 
or O yz quadrupoles, this would be one of the components of J, which should be 
easy either to check or to rule out. However, a pre-existing octupolar symmetry 
can also be broken by a transition to a new octupole-hexadecapole order, or by 
a combination of 0\ quadrupole and triakontadipole order. 

It is interesting to notice that if recent NQR results 1 ' on URu2Si2 are inter- 
preted as a hint that the onset of octupolar hidden order at Tq = 17K is followed 
by quadrupolar ordering at T* « 13. 5K, this sequence of events may fit several 
of the scenarios found in our general classification scheme. However, we have 
to await further evidence showing that the NQR anomalies at T* sa 13.5K are 
associated with an equilibrium phase transition. 

§1. Introduction 

There is increasing interest in orbital ordering phenomena, and their relationship 
to magnetism EJ>EJ>EJ Many of these phenomena can be at least partially understood 
in terms of localized electron models. This is clearly justified for Mott-localized d- 
electrons. Some /-electron systems are semiconductors (like NpC^), but more often, 
we wish to describe multipole ordering in metals like rare-earth-filled skutterudites 
or URu2Si2- In the standard modeP* of the majority of rare earth elements the 
/-shell is localized, and we may invoke a similar feature for the systems of interest 
to us. It can be argued that itineracy and multipolar ordering are complementary 
features, and that they may be manifest in different phases of the same /-electron 
system. In this case, the localized description is acceptable for the ordered phases, 
even though we know that it could not be extended to the whole phase diagram. 

The common starting point is the existence of an iV-dimensional local Hilbert 
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space \n), which allows the definition of N 2 local operators \n)(m\ (of these, N 2 — 1 
are non-trivial, and are also called local order parameters). N(N + l)/2 operators 
have symmetrical, and N(N — l)/2 operators have antisymmetrical character. They 
can be chosen as 

X s mn = \{\m){n\ + \n){m\) (1-1) 

for symmetrical and 

X a mn = ^{\m){n\-\n){m\) (1-2) 

for antisymmetrical operators where i was inserted to ensure hermiticity. 
There are two canonical cases: the basis may consist of 
I. N/2 time-reversed pairs (the N-fold degeneracy arises as 2x(N/2), 

(Kramers degeneracy) x (non-Kramers degeneracy) 
II. real basis states each of which is time reversal invariant (the N-fold degeneracy 
is purely non-Kramers degeneracy) 

Case I. The local Hilbert space consists of {N/2) time-reversed pairs. There 
are N(N — l)/2 time-reversal-even order parameters (this includes the trivial 1), and 
N(N + l)/2 time-reversal-odd order parameters. 

N = 4 is realized by the Jg representation of the cubic double group, the irrep 
of the ground state level for either the 4/ 1 compound CeBg (RefsP®), or the 
5/ 3 compound NpC>2 (RefP). Six operators: i and the five quadrupoles are time- 
reversal-even, while ten order parameters: the three dipoles and seven magnetic 
octupoles, change sign under time reversal. A general form of intersite interaction is 
a sum over the quadratic invariants, for a cubic system with six independent coupling 
constants (one dipolar, two quadrupolar, three octupolar). It was, however, soon 
realized^ that the consideration of the SU(4) symmetrical model (with all couplings 
set equal) should be enlightening. Later research showed that CeBg can be regarded 
as a "nearly SU(4)-symmetrical" systemP^JJnJ 

Case II. If all basis states can be chosen real, the X^ n are real operators 
and the X^ ^ are imaginary operators. Choosing one of the real operators as the 
projection onto the entire local Hilbert space 1 = \m){m\, and orthogonalizing 
the remaining N — l diagonal operators, we are left with N(N+l)/2— 1 time-reversal 
invariant local order parameters. N(N — l)/2 local order parameters (the imaginary 
operators X^ n ) change sign under time reversal. 

This case is realized for non-Kramers ions which have an even number of elec- 
trons. We are interested in the 5/ 2 (U 4+ ), and 4/ 2 (Pr 3+ ) configurations as they 
appear in URu2Si2 and PrFe4Pi2, resp. It is always debatable which N to choose. 
Crystal field doublets (N = 2) have been suggested for both systems but it is unlikely 
for PrFe4Pi2, and at least not widely accepted for URu2Si2- In fact, for both systems 
we settle for a local Hilbert space which is composed of the bases of several irreps. 
For PrFe4Pi2 the quasi-quarte r± + i~4 seems to work. Under tetrahedral sym- 
metry, the nine time-reversal-even order parameters are a singlet, five quadrupoles, 
and three hexadecapoles, while the six time-reversal-odd order parameters are three 
dipoles and three octupolesP^ 
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1.1. The nature of the higher multipoles 

While dipoles and quadrupoles are well-known from decades of research experi- 
ence, it is only recently that octupoles were seriously considered, and higher multi- 
poles are virtually never discussed. It is worth forming a picturial idea of octopule-, 
hexadecapole-, and triakontadipole-carrying shells. 

The J = 4 Hund's rule two-electron states of U 4+ = 5f 2 ionic cores are rather 
complicated objects, superposed of a number of IjSj components by the Clebsch- 
Gordan coefficients. However, there is strong similarity between the multipoles com- 
posed of two /-orbitals and two spins via a number of projections, and simplified 
multipoles built of I = 4 orbitals (atomic ^-states). We derive the charge cloud 
shapes and current distribution patterns from acting with the analogous purely or- 
bital (e.g., J Z (J 2 ~ Jy) ~^ lz(lx ~ ly)) operators on fictitious atomic states (for our 
purposes, fourth-order spherical harmonics 4>) for ^-states). In this Section, 

we do not discuss crystal field effects. 

We diagonalize multipole operators in the 9-dimensional / = 4 space, and base 
our pictures on the state with the (in absolute value) largest eigenvalue. This can 
be interpreted as finding the ground state of an effective field which orders the 
multipole. We mention that free-ion multipole eigenstates are not necessarily the 
same that we find within the restricted Hilbert spaces of crystal field problems. This 
will be illustrated by comparing Fig. |2] and Fig. [5J Furthermore, smaller Hilbert 
spaces may not give a representation of certain multipoles at all. 

For odd-rank (magnetic) multipoles, the eigenstates appear in time-reversed 
pairs, similar to the l z = ±m dipole eigenstates, and the multipole spectrum is 
symmetrical about 0. For even-rank (electric) multipoles, each eigenstate can be 
chosen time-reversal invariant. 

1.1.1. Octupoles 

First, let us consider the octupole 

1 —i 

T~xyz * ~^lxlylz = (l-\-lzl— l—lgl—J ( 1 '3) 

The simplest octupolar states are realized in the I = 2 subspace of e g functions 

*±(Z = 2) = ^= (\3z 2 - r 2 ) ± i\x 2 - y 2 )) (1-4) 

<&±(l = 2) form a time-reversed pair. 

The charge and the current distribution of ^+(Z = 2) is shown in Fig. ^ The 
object got its name from the eight magnetic poles: of the eight current eddies we see 
in the figure*), four belong to magnetic field lines entering, and four to those leaving 

*^ Two general remarks about Figures 00 charge density angular dependences are shown, 
i.e., the lobe shapes do not include the radial fall-off of the atomic wave functions. Flow lines are 
to indicate the sense of circulation of the current but (in contrast to the textbook interpretation) 
the density of flow lines is not associated with higher current density, but is rather arbitrary. Flow 
lines were calculated by solving the differential equation for the tangential curves for the calculated 
vector field of currents, and initial values were randomly generated. (We found that the direct 
representation of the vectors would give unattractive figures). 
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Fig. 1. Left: Charge distribution (lobes) and current pattern (flow lines) in the T xyz octupole 
eigenstate in the I = 2 Hilbert space. Right: the normal component of the magnetic field is 
pointing inwards (outwards) in the light grey (dark grey) areas. 

the surface. Note that the magnetic field pattern measured^ by /iSR in NpC>2 bears 
an overall similarity to what we expect in the neighborhood of an octupole-moment 
bearing shell. For the time-reversed partner = 2) we would find the same charge 
cloud with reversed currents. 

The concept of octupolar ordering was pioneered by Korovin and Kudinov^ 
who envisaged an antiferro-octupole pattern of the e g states p-4|) resulting from spin- 
orbital exchange in Mott insulators. A similar possibility arises in the E doublets 
of trigonal compounds^ where, however, the octupole moment is mixed with the 
orbital moment l z . In any case, antiferro-orbital order is likely to be combined with 
spin ferromagnetism!^'^'^ Metallic phases, including the itinerant octupolar 
phase, were investigated for the e g band Hubbard model by Takahashi and Shiba.^ 

Currently known realizations of octupolar order appear in /-electron systems. 
Field-induced octupoles play a role in understanding the phase diagram of CeE>6 
(RefP). Kusunose and Kuramotd^ called attention to the fact that octupole mo- 
ments are ideally suited for the role of "hidden" primary order parameters. A de- 
tailed study by Kubo and Kuramoto^ makes a convincing case that the "Phase 
IV" of Cei-^La^Bg is an antiferro-octupolar phase. At about the same, it became 
accepted that the long-standing mystery of the nature of the 25K transition of NpC>2 
is solved by identifying it with the triple-q ordering of i~5 U octupolesPSJ^SJ''^ A 
symmetry analysis was successful in identifying the unique octupolar signature in 
NMR spectra on NpC-2, and Cei_ x La x Be (RefP^ The phenomenology of the two 
systems shows certain similarities, as it is to be expected since, e.g., the relationship 
between the anomalies of the linear, and the third-order, susceptibilities follow from 
general thermodynamic reasoningP^''^ The same should be true of URu2Si2 whos 
hidden order, we argued, 2,1 is also of octupolar nature. 

Fig. ^ illustrates the symmetry of the (uniform) octupolar ground state with 
{T X yz) as the order parameter. We leave to Sec. I3,2l the construction of an octupolar 
symmetry group, which will be carried out for a combination of tetragonal crystal 
field and octupolar effective field, the case relevant for URu2Si2- Here we use Fig. ^ 
to visualize the hybrid nature of some of the symmetry elements. The charge dis- 
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tribution is highly symmetrical (octahedral). However, the true symmetry is that 
of the magnetic field pattern, so in purely geometrical terms it is tetrahedral. It 
is obvious, however, that C4 rotations are allowed if they are combined with time 
reversal T (which reverses the currents and the field). Including elements like C4T, 
a higher non-unitary symmetry group could be derived. 

Figure. H shows why octupolar order can be so well "hidden". The charge dis- 
tribution (which influences directional bonding, i.e., the structure of the crystal) is 
cubic, hiding the fact that if current directions are considered, the pure C4 axes are 
not symmetry elements. Similar considerations arose when it was asked whether the 
pseudocubic phase of Sr-doped LaMnC>3 is not, in fact, octupolar ordered.^ 

One might have thought that the character of the octupoles is the same whatever 
Hilbert space we represent them on, but there are interesting nuances here. 




Fig. 2. Current distribution in the T xyz octupole ground state (|1 ■ 5|> in the 1 = 4 Hilbert space. 
Only one octant is shown. 

Within the I = 4 shell (which mimicks*) URu2Si2's 5/ 2 — > J = 4 multiplet), 
T xyz has the minimal eigenvalue — 6v3, corresponding to the eigenstate (expressed 
in the \l z ) basis) 

^ = ^(1 + 4) + I - 4)) - 2iyJL(| + 2) + I - 2)) - yi|0) (1-5) 

(the maximal eigenvalue 6\/3 corresponds to the time-reversed of ^Q Ct ). The current 
distribution for <pQ Ct is shown in Fig. |2 The overall pattern of eight alternating 
vortices is the same as in Fig. ^ However, within each of the major eddies four new 
sub-eddies appear: three (arranged like the petals of a flower) rotating in the sense of 
the eddy as a whole, and a little central eddy counter-rotating. The symmetry of the 
current distribution in Fig. |2]is the same as that in Fig. ^ but there are differences 
in the magnetic field pattern. While in Fig. ^ the magnetic field is maximum in 
the center of an octant, and has the same sign everywhere within the octant, in 
the I = 4 solution (Fig. I2J the field changes sign in a small central region of the 

*' However, a calculation of the currents for realistic two-electron states would be desirable. 
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octant and appears smaller than in the "petals". On the whole, the magnetic field 
of the octupolar currents is weaker for I = 4 states (f 2 configurations) than for the 
simplest I = 2 solution. This may account for the weakness of the internal fields in 
URu2Si2- It should be noted, though, that the current distribution shown in Fig. [2 
was obtained for a free ion subject to the octupolar effective field only. The situation 
changes if crystal field effects are included (Sec. 13. 2j) . 

It would be of obvious interest to calculate internal field distributions for oc- 
tupoles, and maybe also for triakontadipoles, for these predictions could be checked 
against /xSR and NMR observations. We note that Kubo and Kuramoto 211 discussed 
the situation for the two well-established octupolar systems NpC>2 and Cei_ x La z B6. 
for which 500G and 40G, respectively, are measured; both are in excess of the theo- 
retical estimate. The measured internal field in URu2Si2 is much weaker: 29 Si-NMR 
linewidth gives < 12G_(RefP), while at the muon stopping sites, the internal field 
is at most 1-2G (RefP). 

1.1.2. Hexadecapoles 

The hexadecapole l x ly(lx — ly) has the largest eigenvalue for the eigenstate 

aW = -^(|2)-i|-2>). (1-6) 

In spite of the appearance of i, ^hex is time reversal invariant. It can be chosen real 
and its positive and negative lobes can be identified (Fig. |3J). There is a general 
similarity to quadrupolar eigenstates (which also have positive and negative lobes), 
only the number of lobes is higher. We do not construct the symmetry group of the 
hexadecapole effective field but we notice that, like the octupolar symmetry group, 
it would have composite symmetry elements: purely geometrical transformations 
combined with "lobe sign reversal". 

Hexadecapoles as order parameters were discussed within the tetrahedral sym- 
metry classification,^ valid for Pr-filled skutterudites. We are not aware of the 
existence of primary hexadecapole order in any system. 

1.1.3. Triakontadipoles 

There are 11 triakontadipoles but we consider only l x l y l z (l 2 — ly) which appears 
as the lowest-rank A\ u multipole in the tetragonal classification (TableHJ). Its ground 
state is 

£tria = ^(|4) + |-4)+n/2|0)) . (1-7) 

with eigenvalue 6\/35. As shown in Fig. El there are 32 cells of alternatingly flowing 
current, so the associated magnetic field is more short-ranged than in the case of 
an octupole (Fig. As in the case of octupoles, the symmetry of the charge 
distribution is higher than that of the current distribution: the charge cloud has a 
Cs axis which is reduced to C4 for the currents. However, the combination of the 7r/4 
rotation with time reversal C%T is a symmetry operation. It would be interesting to 
meet triakontadipolar order in nature. 
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Fig. 3. Charge distribution in the electric hexadecapole eigenstate 11-61 . The sign of the lobes of 
the wave function is also shown. 




Fig. 4. The charge (left) and the current (right) distribution in the magnetic triakontadipole eigen- 
state im . 

§2. Recent experimental developments on URu2Si2 

The magnetic behavior and phase diagram of the intermetallic compound URu2Si2 
have been the focus of attention for over two decades Specific heat measure- 
ments' 50 ' show that electronic entropy of 0(ln2)-/cB is released by the time the tem- 
perature reaches 30K, and a sizeable fraction of it is associated with the A-anomaly 
at To ~ 17K. Thus it is not unjustified to think of the phase transition as the full- 
scale ordering of a localized degree of freedom, but the nature of the order parameter 
remains hidden. It is obviously not the tiny (~ 0.03/iB) antiferromagnetic moment 
which is observed by neutron scattering.^ 

Small (either static or slowly fluctuating) moments have long been held to be an 
attribute of heavy fermion systems on the borderline between localized and itinerant 
/-electron phases. The specific heat value would allow to classify URu2Si2 as a "light 
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heavy fermion system", raising the question whether its mysterious properties may 
be related to an exotic itinerant phase of strongly correlated /-electrons. Orbital 
magnetism due to plaquette currents,^ unconventional density wavesJ^J^ and 
Pomeranchuk instability leading to a nematic stated are in this category. 

The RKKY interaction mediates a variety of multipole-multipole interactions 
between the /-shellsP^ However, the ordering may be foiled by the collective Kondo 
effect: the formation of a heavy Fermi sea with a large (Luttinger) Fermi surface 1 ^ 
may swallow up the localized moments. The Kondo-to-RKKY transition has been 
extensively studied for the case when the relevant local degrees of freedom are spin 
dipoles 38 ' but work on the general multipolar problem has only just started.^ It 
is an intriguing possibility that the 17K phase transition in URu2Si2 coincides with 
the itinerant-to- localized transition of the /-electrons!^ This would lend credence 
to describing the T < Tq order in terms of a localized /-electron model even if a 
more satisfactory description will have to encompass itinerant aspectsPJ 

It has long been knowrPS'EU that the symmetry classification of local order 
parameters in tetragonal crystal fields follows Table [fl For each of the irreps, only 
the lowest-rank multipole is listed. The result can be viewed as arising from the 
tetragonal splitting of the operator irreps of the cubic classification schemed 

Table I. Symmetry classification of the local order parameters for B = (X>4h notations^ overline 
means symmetrization) . g (u) belongs to even (odd) parity under time reversal. 



sym (.9) 


operator 


sym (u) 


operator 




£ 


A lu 


JxJyJzi^Jx Jy) 


Big 


1~L\ — JxJy^Jx Jy) 

ol 


A 2u 
Biu 


J, 

Txyz — JxJyJz 


B 2 g 
Eg 


&xy — JxJy 
{&xz •> &yz } 


B 2u 
E u 


Tz ~ Jz(J& ~ Jy) 
{Jx, Jy} 



Progress with uncovering the true nature of the hidden order of URu2Si2 has 
been particularly slow because of a seemingly extrinsic property of almost all of 
the samples: they show a kind of micro-antiferromagnetism with moments m ~ 
0.02 — 0.04/ie || z with the simple alternating order Q = (0, 0, However, 
evidence from microscopic measurements^ EDEIJIEJ shows that the nominally small 
moment should be understood as a relatively large moment at a minority of the sites. 

The relationship between the hidden order parameter ifi and the antiferromag- 
netic moment m has been, and still remains,^ a matter of debate. There are two 
basic possibilities A) hidden order and antiferromagnetism share the same sym- 
metry,^ or B) they are of different symmetries. In Case A), the relative amplitude 
of m and ip can be tuned continuously (e.g., by pressure), and there is no sharp 
distinction between the two orders. In Case B), the two orders are incompatible, 
and there must be a first order transition from the ip 7^ 0, m = phase to the 
ifj = 0, m ^ phase. In the present discussion (as in RefP) we take the conclusion 
drawn from recent /^SR experiments^ as our starting point: there is a first-order 
transition, the symmetry of ip is different from that of m, and therefore the hidden 
order parameter must be sought from among other entries in Table [|J 
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As far as gross features like linear and non-lineaJ^ susceptibility, specific heat, 
etc. are concerned, with a little adjustment of the crystal field level scheme, B\ g 
and F>2 g quadrupolar, 52 ' and B\ u and B2 U octupolar^' models do about equally 
well. Evidence beyond this simple range of experiments has to be invoked to choose 
between the quadrupolar and octupolar scenarios. 

We cite two crucial (but as yet unpublished) experiments to argue that the 
hidden order is not quadrupolar. First, there is a remarkable mechanical-magnetic 
cross-effect (at T < 17K): in the presence of uniaxial stress applied perpendicularly 
to the tetragonal main axis, large- moment antiferromagnetism with the simple pat- 
tern described above becomes visible for neutron magnetic Bragg scattering. 53 ' This 
is unambiguous proof that the background order breaks time reversal invariance, and 
is thus certainly not quadrupolar; the simplest remaining choice is octupolar order. 

The second evidence is coming from recent NQR measurements by Saitoh et alP 
The temperature dependence of the electric field gradients was followed carefully 
from relatively high temperatures (70K) to well below To = 17K. V zz changes all 
the time, reflecting the T-dependence of the tetragonal crystal field component 0\. 
Saitoh et al's findings can be formulated in two statements: (a) Neither of the field 
gradients shows anomalies at To = 17K, the onset temperature of hidden order, 
so the hidden order cannot be quadrupolar. (b) There is an anomaly in the NQR 
signal at T* as 13. 5K, so something happens to the quadrupoles there. A possibility 
is that T* is the ordering transition of quadrupoles; then it has to be a second HO 
transition following the first one at To. We emphasize that a lot more experimental 
evidence (especially specific heat, susceptibility, etc) is needed before the existence 
of a phase transition at T* may become accepted. It is not our aim to describe 
the two transitions in any detail. We merely emphasize that, given the symmetry 
of URu2Si2, octupolar ordering can be followed by quadrupolar ordering; but then 
further induced magnetic multipoles should be observable. 

While circumstantial evidence for octupolar ordering looks encouraging, at- 
tempts for its direct verification have as yet yielded negative results. An early, 
very specific neutron scattering investigation ruled out either B\ u or B2 U octupoles 
at the selected wavevectors, including the (0,0,1) periodicity of the weak-moment 
antiferromagnetism^^ It has been argued that resonant X-ray scattering would not 
see the octupoles.^ Last, but not least, though the octupolar scenario is as yet 
the only one to account for the observed fact that transverse uniaxial stress induces 
|| z antiferromagnetism, in a simple mean field versiorP it forces a choice between 
a || (100) and a || (110) directions of stress, while experiments seem to tell us that 
these directions are equivalent.^ An essential insight is missing here. 

A microscopic theory will have to address the questions raised above. Our 
present investigation is of a limited scope: we use general arguments to classify the 
symmetry- allowed ordering transitions of URu2Si2- We are particularly interested 
in the possibility of a sequence of such continuous phase transitions. We have no 
statement to make if the lower-temperature transition is of first order. 

We address the simplest questions for which no knowledge of microscopic details 
is required: Once we had an octupolar ordering transition, is there any compelling 
reason to expect a second symmetry-breaking transition? What may be the order 
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parameters? Should we expect still more, as yet undiscovered, phase transitions? 

In an abstract sense, our question is the following: assuming the symmetry of the 
high-temperature ("para") symmetry group G pa , Ta has been broken by a spontaneous 
ordering transition, introducing order {0\) ^ lowers the symmetry to Q(0\) £ 
<5para- What is the structure of Q{0\)1 What is the new classification of the order 
parameters? Assuming that such order parameters are found, does this imply that 
further continuous phase transitions necessarily happen? 

In Ref.P we discussed in some detail two questions of this nature: the symmetry 
classification of the order parameters in the presence of an 1) external magnetic field 
B || z, and 2) a uniaxial stress. Formally the same question arises if instead of 
externally applied fields, we assume the presence of an effective field associated with 
either a dipole ordering transition (1) or quadrupolar order (2). 

In Sec. 13. 11 we rephrase our earlier results on the effect of an external magnetic 
field, and add some remarks about the field direction dependence. In Sec. 13.21 we 
analyze the symmetry in the presence of a T xyz octupolar (effective) field, and the 
possibility of a sequence of phase transitions. Sec.0]illustrates the general arguments 
with a simple mean field calculation. 

§3. Symmetry lowering in external, or effective, fields 

3.1. Magnetic field 

In the presence of an external magnetic field H z || (0,0,1) or equivalently 
(J z ) > 0, the remaining purely geometrical symmetry elements are £, 2C$, and C\ 
(and, naturally, the inversion I). The geometrical symmetry (described by unitary 
operations) is lowered to C^h, but the full*-* symmetry group Qteti(Jz) contains also 
non-unitary elements, namely 2C' 2 T, and IC'^T , where T is time reversal. Qtetr(Jz) 
contains eight elements, its character table (and multiplication table) is the same as 
that of the tetragonal point group P4. The character table is given in Table ITTl 

Reducing the symmetry from Q te tr to Qtetr (Jz), some of the originally inequiv- 
alent irreps of Qtetr become equivalent. This parentage of the irreps is shown in 
Table [HI It also follows that the corresponding order parameters get mixed. This 
was discussed in RefP 

Table II. Character table of the non-unitary symmetry group <5tetr( Jz) of a tetragonal system placed 
in a magnetic field B = (0, 0, B z )- Space inversion is omitted. 
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Ai 
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1 


1 
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1, Jz 


A 2 
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-f 


-f 


Alu, A 2 g 


JxJy^Jx Jy^i JxJyJz{^Jx 


Bi 


1 


-1 
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1 


-f 


Big, B 2u 


OlTf = Jz{Jl~Jl) 


B 2 


1 


-f 


f 


-f 


f 


B 2 g, Blu 


Oxyi Txyz 


E 


2 





-2 








Eg, E u 


{ Jx, Jy } , { Oxz, Oyz } 



*' Magnetic field is invariant under space inversion. It is understood that all symmetry groups 
would get doubled if we included the inversion T. 
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a; is a lower-symmetry direction than z and correspondingly a system subjected 
to a magnetic field B = (B x , 0, 0) has a smaller symmetry group <5tetr(A0 ( Table ITTTf) . 
In contrast to the case B || z where irreps taken from the zero field case kept 
their dimensionality only pairwise merged, now the two-dimensional irreps split. 
The parentage of the B || x irreps is as follows Ai g ,Bi g — > A, A2g,E>2 g — ► A, 
Ai u , Biu — ► A, A2u, B2u — > A, E g — ► A + A, E u — ► A + A- 

Table III. Character table of the non-unitary symmetry group Gtcti{Jx) 



irrep 
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CfT 


C*2x 




A 


1 


1 


1 


1 


r 2 


1 


1 


-1 


-1 


r 3 


1 


-1 


1 


-1 


A 


1 


-1 


-1 


1 



The presence of non-identity irreps in Tables [H] and II I II shows that for fields 
B || z, and even for B || x, a variety of symmetry breaking transitions are possible. 
Since the field has broken time reversal invariance, and induced a polarization in its 
direction*^ , all these transitions have the character of orbital ordering (lifting some 
residual non-Kramers degeneracy), and also lead to the appearance of time-reversal- 
odd transverse polarization components. We enumerate the possibilities (B || z) 

• A field B || z allows the ordering of E type quadrupoles and this is accompanied 
by the appearance of dipole polarization in the xy plane. This was suggested^ 
for the disjoint high-field phase observed in experiments^ on URu2Si2- Anal- 
ogous phenomena are observed 56 ' in high field experiments on the tetrahedral 
skutterudite PrOs4Sbi2- 

• The ordering is of B\ type. In our suggestion^ for the low-field order of 
URu2Si2, the primary order was Fz octupolar (this survives when the field 
is switched off), and 0\ quadrupoles were field- induced. Santini's scenario^ 
was the opposite: primary quadrupolar order, and field-induced octupoles. 

• B2 order parameters: primary F xyz octupolar order and field-induced O xy 
quadrupoles. In the present paper, we assume this is the hidden order with 
onset temperature Fq = 17K. Lacking a microscopic model, it is impossible to 
decide between this scenario and the previous B\ scheme. 

• The exotic possibility of A2 order: TL\ hexadecapoles and J x J y J z {Jx ~ Jy) 
triakontadipoles, has not been considered yet. 

The discussion of the cases B || y and B || (x ± y) would be analogous to 
B || x: the remaining symmetry group has order 4 (or 8, if we include inversion). 
The symmetry group for fields lying in the x~y plane in non-special directions is 
generated by C\F and F. For fields in general out-of-plane directions, the only 
remaining symmetry is inversion. 

We conclude that a spontaneous symmetry breaking transition in external fields 
remains possible if the field is either parallel or perpendicular to z. In particular, the 

*' Assuming the local Hilbert space allows a Zeeman splitting for the particular field direction. 
E.g., the quasidoublet {jti), ^4}} of URU2S12 would not split in a field B || z. 
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octupolar transition can remain second order for a range of field strengths. However, 
T xyz octupoles get mixed with O xy (B2) quadrupoles for B || z, and with ^-derived 
quadrupoles (a suitable combination of O yz and O xz ) for Biz. For general field 
directions, a smearing of the octupolar transition should be observed. This effect 
may be helpful in deciding whether the hidden order is of octupolar nature. 

The above reasoning can also be extended to discuss further symmetry breakings 
after a phase transition to a dipole-ordered phase has taken place. The effective field 
of the ordered moments acts the same way as an external field. E.g., we may conclude 
that following a transition to a (J z ) ^ phase at the higher critical temperature T> , 
tranverse dipoles (from the E doublet { J x , J y }) may also order at a lower critical 
temperature Tf-, and that (J x ) ^ induces (O xz ) ^ 0. Though symmetry allows 
this (or some other scenario identifiable from Table [nj to happen, whether this 
potentiality is realized depends on the nature of the microscopic model. 

3.2. The symmetry group of the octupolar phase 

We may ask whether following the onset of, say, T xyz type octupolar order a 
further symmetry lowering transition is possible. This is the same question as to 
whether spontaneous symmetry breaking in an external octupolar field of B\ u sym- 
metry is possible. The question is partly academical, but it is also motivated by 
the NQR findings on URu2Si2 by Saitoh et alP In particular, they claimed that 
either O xy or 0\ quadrupoles appear at T* « 13. 5K, well below To « 17K which we 
associate with the onset of octupolar order. However, as we emphasized before, the 
connection of our arguments with experiments remains tenuous. 

We have to identify which operations leave the T xyz octupoles unchanged, i.e., 
we are looking for the octupolar symmetry group Gtetr{T xyz ) a s a subgroup of Gtetr- 
Schematically, T xyz is represented in Fig. |S] where the current distribution in the 
ground state of T xyz is shown for two- and three-dimensional subspaces selected 
from the 9-dimensional Hilbert space of the tetragonal crystal field eigenstatesP**^ 
It is obvious that the selection of the basis functions (which is "done" by the crystal 
field potential) has a strong influence on the details of the current pattern. However, 
though the pattern is much more decorated for the ^3)} doublet than for the 
triplet { 1 1\), ^2)) 1*4)} (the crystal field model used in RefP), the symmetries are 
the same. In either case, there is a current distribution with eight major eddies: 
four positive, and four negative (whether within these, there are sub-eddies, has 
no influence on the symmetry). Though there are local magnetic fields, the total 
magnetic moment is zero. A C4 rotation takes positive eddies into negative ones, 
and vice versa; however, reversing also the direction of currents, the original state is 
restored. From the C' 2 and C 2 elements of (/tetn the latter two have to be combined 
with T . The character table of the symmetry group Gtetr(%cyz) is shown in Table llVl 
(as before, it is understood that the space inversion T would generate the other 
half of the complete symmetry group). We have also given the resulting symmetry 
classification of some of the order parameters in the last column of Table IIVI 

It is interesting to note the similarities and dissimilarities to the symmetry group 
of J z f Table Hi)) . The character table is the same only in Gteti(%cyz), 2C4 and 2C 2 , 
while in Gtetr(Jz), 2C 2 and 2C 2 have to be combined with T. 
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Fig. 5. Crystal field effects on the current distribution in the T xyz octupole ground state. Left: the 
doublet {», |t 3 » (|ts) = (l/\/2)(|2) + | - 2))). Right: the triplet {\ti}, \t 2 ), \U}} 
The two current distributions have the same symmetry given in Table llVl 

Table IV. Character table of the non-unitary symmetry group Gtetr(1xyz) of a tetragonal system 
placed in an effective octupolar field. Space inversion is omitted. 



irrep 


£ 


2C 4 T 


cl 


2C' 2 


2C^T 


parentage 


operators 


A x 


1 


1 


1 


1 


1 


Ai 9 , B lu 




1 T 


A 2 


1 


1 


1 


-1 


-1 


A 2 g, B 2 u 


JxJyi, Jx 


J%f)i Tz ~ Jzi^Jx Jy) 


Bx 


1 


-1 


1 


1 


-1 


Big, Ai u 




JxJyJzi^Jx Jy ) 


B 2 


1 


-1 


1 


-1 


1 


B 2 g, A 2u 




Qxyi Jz 


E 


2 





-2 








Eg, E U 




Jy } > { Oxz, @yz } 



Analogous results would have been obtained if we had assumed a T" octupolar 
field instead of T xyz . That was our starting assumption in RefP 
GtetriXcyz) has two generators. There are several choices: 

• C{T and C' 2x or C' 2y 

• C 4 T and C'l x+y T or C' 2 '^ y T 

• either of C' 2 and either of C 2 

At T* < To, one of the order parameters appearing in Table ITVl acquires non- 
zero expectation value, and the symmetry of the system is lowered to one of the 
subgroups of Gtetr{%cyz)- The list of the subgroups is 



Qi 


— {£ ,C/iT ,Cl,C^T} 


G-2 


= {£AJ 


03 


= {£A,y} 




= ^2,x+y^~} 








= {^i^2,xi C%yiC\} 


Qi 


= {^iC 2x+y T, C 2x _ y T, Cf} 


Gs 


= {£,cj} 



Each of the order parameters appearing in Table llVl breaks one, or several, of the 
symmetries in Gtetr {Xcyz), and thus reduces the symmetry to one of the subgroups. 
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All the possibilities are listed below 



T/ and J x J y {Jl - J2) — & 



and Jr — ► £ 2 
C X2 and J y — ► G 3 (3-1) 



2 2 and J x JyJ z (Ji-J2)^G 6 

O xy and J z — ► £7 (3-2) 



Time-reversal-even and time-reversal-odd order parameters appear in pairs. Since 
the octupolar background already breaks time reversal invariance, it cannot be "un- 
broken" by the next transition, so the new phase has a new pattern of non-zero 
current densities. All the quadrupoles allowed by tetragonal symmetry can appear 
in a symmetry-breaking transition but they bring either dipoles or higher magnetic 
multipoles with them. We list the possibilities: 

• Most straightforward is the case of O xy quadrupolar ordering accompanied by 
magnetic moments in the z direction. The coupling of these three multipoles: 
T X yz, Oxy, an d J z has been discussed from several points of view. In Ref.P 
we argued that on the background of T xyz octupole ordering, applying uniaxial 
stress a || (1,1,0) creates O xy quadrupoles, and therefore magnetic moments J z , 
offering an interpretation of the experimental results by Yokoyama et al. 53 ' The 
existence of a third-order invariant T xyz O xy J z implies that if in an ordered phase 
T xyz 7^ then also the correlator (O xy J z ) 7^ 0. The non- vanishing correlator 
does not automatically give non-zero values to (O xy ), and (Jz), but at least it 
gives a hint that O xy , and J z are likely actors in a cooperative phenomenon. 
Here we discuss them as coupled order parameters of a low-temperature phase 
transition following a high-temperature octupolar transition. 

• O xz (O yz ) quadrupoles would be accompanied by the transverse magnetization 
component J x (J y ), so the experimental verification (or refusal) of this scenario 
should be straightforward. 

• More exotic is the possibility of 0\ quadrupolar ordering (one of the likely can- 
didates according to RefP) which should be accompanied by J x J y J z (J% — J^) 
triakontadipole ordering. 

• Finally, the symmetry of the T xyz octupolar field can be spontaneously broken 
by ordering the T z octupoles; this should be accompanied by hexadecapole 
order of the Hi = J x J y (J x — J y ) type. We note that TLi hexadecapoles at 
the U sites can combine to an effective quadrupole at the Ru sites, so this 
scenario is not necessarily in conflict with the findings by Saitoh et al. 1 ' It is 
an added attraction that the simultaneous presence of (T xyz ) 7^ and (Tf) 7^ 
would allow that uniaxial press applied either in the a || (100) or the a || (110) 
direction induce J z 7^ 0, as observed.^ We note that the measurements of 
Yokoyama et alP^ were carried out at 1.4K, well below either To or T* , so both 
octupolar amplitudes would be near their saturation values. We have to admit, 
though, that in our scenario it would be difficult to get them equal. 
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3.3. Octupolar phase in external magnetic field 

It is of some interest to combine the previous two cases to discuss the remaining 
possibilities of symmetry breaking if the established T xyz octupolar order is subject 
to an external magnetic field. 

A field B || z reduces the symmetry to the four-element subgroup Gj. There 
are four irreps: A[, A 2 ,A' 3 , and A'^. In terms of the irreps shown in Table ITVl their 
parentage is: 

A l ,B 2 ^A' l , 
A 2 ,B 1 ^A' 2 , 
E^A^ + A' 4 

Thus now T xyz , J z , and O xy are all present in the identity representation. 

A spontaneous symmetry breaking transition is possible to a phase with (Of) 7^ 
0; but then Tf, Hi, and J x J y J z (J^ — J|) are induced order parameters. 

There are two more possibilities of ordering, both involving a transverse dipole 
and a quadrupole (e.g., J x — J y and O xz — O yz ). 

A field B || x reduces the symmetry to the two-element subgroup Q%. Still, there 
remains one symmetry element to break: it can be done by J z , and a number of 
associated multipoles. 

To conclude this subsection: if the T < Tq hidden order is octupolar, the pos- 
sibilities of remaining symmetry breaking at T* depend very much on the direction 
of the applied field. It seems that the only direction offering non-trivial possibilities 
is B || z, where the 0\ quadrupolar symmetry breaking is accompanied by induced 
octupolar, hexadecapole, and triakontadipole moments. 

§4. Mean field calculations 

From the fact that the symmetry group of a model contains non-trivial elements, 
it does not necessarily follow that symmetry breaking phase transitions occur until 
the only remaining symmetry element is the unit operator. It is possible that the 
local Hilbert space is not large enough, or it does not have the right structure, to 
support a sequence of ordering transitions. An ordering transition is an inevitability 
only if in its absence, the ground state is degenerate. In models of octupolar ordering, 
this is usually not the case. Santini's^S and ouJ^ crystal field model is built with 
three singlets, so the single ion ground state is in any case non-degenerate. Still, 
symmetry breaking transitions can occur by the induced moment mechanism if the 
level splittings are not too large to begin with. After the first ordering transition - 
which we assume is T xyz octupolar ordering - has taken place at To, the symmetry is 
reduced to Gtetr (%cyz)- The ionic level scheme is formed by the combined action of 
the crystal field potential and the octupolar effective field, and the ground state is 
again non-degenerate. Nevertheless, the induced moment mechanism can be effective 
again and we may ask whether another second order transition may take place at 
T* < Tq, corresponding to one of the options listed in ()3-l|) . A mean field calculation 
shows that this is indeed possible for a wide range of parameters. Naturally, one has 



16 



to assume non-zero couplings for the multipoles which appear in the definition of 
the new order (for instance, 0\ and J x J y J z (J x — Jy), or O xy and J z ), but this is in 
any case more plausible than setting the said couplings to zero. The starting values 
of the crystal field splittings do not play any particular role, except that they have 
to be small enough to allow ordering. 

Let us observe that Gteti(Xxyz) has a two-dimensional (E) representation. It 
shows that if the crystal field scheme includes a low-lying doublet, then this feature 
may be preserved also in the octupole-ordered phase. The two-fold degeneracy may 
then be lifted at the second ordering transition. In this sense, crystal field schemes 
with either a doublet ground state, or a low-lying doublet, offer a more direct route 
to a second symmetry breaking transition. We note that the crystal field level 
scheme is not undisputed: while models with low-lying singlets have been most 
widely discussed, doublet ground state was also considered.^' However, a sequence 
of two continuous transitions is possible either with, or without, a low-lying doublet. 

4.1. Three singlets 

Since f 2 is a non-Kramers configuration, it is possible to choose time reversal 
invariant basis functions. For the singlets A\, A 2 , and I? 2 , these are 

|t 1 ) = 6(|4) + |-4)) + V / l-262|0) (4-1) 

|t 2 ) = ± (|4) - I - 4)) (4-2) 

|t 4 ) = -L (|2) - I - 2)) (4-3) 

In this basis, electric (magnetic) multipoles are real (imaginary) operators (see Ap- 
pendix). 

Choosing the subspace of the singlets (|4-1H4-3|) is the common starting point 
for a mean field calculation in Santini and our** work. The presence of these 
three states is useful for getting good overall agreement with measured macroscopic 
quantities. Depending on the assumptions about the sequence and the splitting 
of the levels, and about the matrix elements connecting them, the three-singlet 
scheme may support J z dipolar, 0\ or O xy quadrupolar, T z or T xyz octupolar, or 
Til hexadecapole order in the first ordered phase, and a combination of two more 
orders from the same list below a second critical temperature. Up to now, mainly 
scenarios with a single transition were considered. SantinP chose quadrupolar order 
while we preferred octupolar order. There is no effective way to choose between Tj? 
and T xyz . Our earlier discussiorP was based on postulating T% ordering. Here we 
phrase our arguments on the alternative assumption that T xyz order appears first. 

The energy (free energy) density can be expressed as a sum of invariants. A 
number of third order invariants contain T xyz . The corresponding third-order con- 
tribution to the free energy is 

AF3 = ci(Qi, Q 2 , Q 3 )WQiWQ2)«i(Q3) 
+c 2 (Qi, Q 2 , Q3)T xyz (Qi)O xy (Q 2 )J z {Q 3 ) 
+c 3 (Qi, Q2, Q 3 )WQi) [MQ2)O gx (Q 3 ) + J y (Q2)O yz (Q 3 )} 



17 



+c 4 (Qi, Q 2 , Q 3 )-WQi)0 2 2 (Q 2 )/C(Q 3 ) (4-4) 

where Qi + Q 2 + Q3 = 0, and K = J x J y J z (Jl - J|). 

If T xyz (Qi) is the hidden order, i.e., {T xyz (Qi)) / at T < To, then in the 
same temperature range the double correlators read off from above become non- 
zero, e.g., {O xy (Q2) Jz(Q,3)) 7^ 0. This is helpful, but in itself not sufficient, for 
the ordering of O xy and J z separately. For that, dipole-dipole and quadrupole- 
quadrupole interactions have to be non-zero, and sufficiently strong to overcome the 
level splittings arising from the combined effect of the crystal field and the octupolar 
effective field. Playing with parameters introduces a lot of arbitrariness into a crystal 
field model, and at present it would be pointless to try to fit them to experiments, 
especially as there is no agreement about the relevant Q vectors. We merely wish to 
demonstrate the possibility of a second ordering transition at T* < Tq. 




Fig. 6. Mean field solution for two consecutive phase transitions. Following the onset of staggered 
T xyz octupolar order, the coupled quadrupolar O xy (Q2) and dipolar J z (Qa) order parameters 
appear at a lower critical temperature (the parameters are given in the text. Note the break in 
the slope of {T xyz (Qi)) at T = T* . This anomaly will be clearly seen in Fig. |7] (right). 

For the order parameters appearing at T* , any of the pairs in Q3-2JI could be 
chosen. For the moment, we arbitrarily pick O xy and J z . 

Ordering is possible by the induced moment mechanism in spite of crystal field 
splittings, but its nature is basically the same as it would be for three degenerate 
singlets*^. We introduced three mean field couplings: X(T xyz ), X(O xy ), and X(J Z ), 
and solved the self-consistency equations for the three order parameters. 

Typical results are shown in Fig. HO The onset of T xyz octupolar order at 
To = 0.645 is followed by the ordering of O xy and J z at T* rj 0.495. To bring 
out the contrast**-* between the complementary parameters O xy and J z , we used 

*^ Crystal field splittings establish an asymmetry between the complementary variables O xy 
and J z because O xy has a matrix element between and | £2} , while J z between \t\) and \t&). 
Asymmetry arises also from unequal multipole couplings. 

**' The critical temperature is the same, but the amplitudes unequal. 
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k„ IV/.(T ) . - . . ■ 

B v xyz B v xyz 

Fig. 7. Left: The temperature dependence of the inverse staggered quadrupolar susceptibility. 
Right: the discontinuity in the temperature derivative of the octupolar moment at T*. 

the parameter set A21 = 1 as the energy unit, A41/A21 = 2, X(O xy ) = A(J Z ), 
X(J z )/A 2 i = 10, and \{% yz )/A n = 17. 

The second step of ordering is assisted by the pre-existing octupolar order. This 
can be seen in the susceptibility plot Fig. (left). We have chosen the quadrupolar 
susceptibility which belongs to the ordering degree of freedom X y{Q2)- The high- 
temperature susceptibility extrapolates to a lower quadrupolar ordering temperature 
than the behavior calculated within the octupolar phase. An alternative argument 
is based on the Landau expansion (we drop the Qs) 

T = a oct (T - T )T 2 + 6 oct T 4 + a quad (T - T x )(0 2 + J 2 ) 

+6 qU ad(0 4 + J 4 ) + eT 2 (0 2 + J 2 ) + fTOJ (4-5) 

The last term corresponds to (|4-4|) . Whatever the sign of /, the sign of OJ can 
be chosen so that the term lowers the free energy, and enhances the lower critical 
temperature from T\ to T*, as observed. 

The octupolar order responds to the onset of quadrupolar-dipolar order by a 
change of the slope of the temperature dependence (Fig. [7J right). 

We emphasize that we did not make any attempt to fine-tune our mean field 
parameters. If we plotted the specific heat, it would show two lambda-anomalies, in 
conflict with the experiments known to us. We think that until further experiments 
are done, it would be pointless to refine our calculation. 

§5. Conclusion 

We studied possible sequences of symmetry breaking transitions in tetragonal f 2 
systems. Assuming that octupolar order sets in first, the next transition may lead 
to mixed dipole-quadrupole, quadrupole-triakontadipole, or octupole-hexadecapole 
order. There is an interesting similarity to the recent NQR results by Saitoh et aP 
but we feel that any detailed comparison would be premature. 
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Appendix 

The off-diagonal multipoles in basis (|4-lH4-3j) are 

-i \ _ / 

J z = 4y/2b ( % 0\ = 2V7 (oOl) (A-1) 

000/ " V 1 



O xy = (-2V146 + 6V5 - 106 2 ) (A-2 
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Hi = 3^35 - 706 2 [lOO 




(A-4) 



All the matrix elements of the triakontadipole operator J x JyJz(J% — Jy) vanish 
in this subspace. 
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